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Abstract Comparative analysis of second harmonic generation in ordinary and 
backward-wave settings is presented. Extraordinary properties of frequency dou¬ 
bling nonlinear optical reflectivity and pulse shaping through phase matching of 
ordinary and backward electromagnetic waves in the nanowaveguides with mixed 
negative/positive spatial dispersion is demonstrated with numerical simulations. 
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1 Introduction 

Metamaterials (MM) are artificially designed and engineered materials, which can 
have properties unattainable in nature. Usually, MM rely on advances in nan¬ 
otechnology to build tiny metallic nanostructures smaller than the wavelength 
of light. These nanostructures modify the electromagnetic (EM) properties of 
MM, sometimes creating seemingly impossible optical effects. Negative-index MM 
(NIMs) are the most intriguing EM materials that support backward EM waves 
(BEMW). Phase velocity and energy flux (group velocity) become contra-directed 
in NIMs. The appearance of BEMW is commonly associated with simultane¬ 
ously negative electric permittivity (e < 0) and magnetic permeability (p < 0) 
at the corresponding frequency and, consequently, with negative refractive index 
n = —y/JH. Counter-intuitive backwardness of EMW is in strict contrast with the 
electrodynamics of ordinary, positive-index materials. Extraordinary properties of 
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nonlinear-optical frequency converting propagation processes has been predicted 
provided that one of the coupled electromagnetic waves falls in the negative-index 
frequency domain and propagates against the others [Tl[51[31[l] . Among them are 
second harmonic generation (SHG) [ll[2j|3l[l] , optical parametric amplification and 
frequency-shifted nonlinear reflectivity [Tl[2l[3l[ll[5] . 


Current mainstream in fabricating bulk NIM slabs relies on engineering of LC 
nanocircuits - plasmonic mesoatoms with phase shifted, negative electromagnetic 
response. Extraordinary coherent, nonlinear optical, frequency converting propa¬ 
gation processes predicted in NIMs have been experimentally realized so far only 
in the microwave [6]. A different paradigm which employs spatial dispersion [71 
[8] was proposed to realize outlined processes. Basic underlaying idea is nanoengi¬ 
neering negative spatial dispersion which would enable negative group velocity 
Vgr = du/dk < 0. Then energy flux, which is directed along the group velocity, 
appears directed against the wave vector. Such a property makes an analog of 
the effective negative refractive index. Various particular realizations of negative 
dispersion were proposed laiioiiniiiiiis]. Basically, many hyperbolic metamate¬ 
rials may support backward-wave electromagnetic modes, see, e.g., na, as well 
as specially designed waveguides |151I16| . This opens new avenues for creation of 
backward-wave photonic devices with unparalleled functional properties. As re¬ 
gards coherent nonlinear-optical propagation processes, critically important is to 
provide for coexistence of ordinary and backward coupled waves which frequen¬ 
cies and wave vectors satisfy to energy and momentum conservation law (phase 
matching). Such possibilities were described in ref. |111I121[T^ . Corresponding MM 
slabs can be viewed as the plasmonic nanowaveguides which support both ordinary 
EMWs with co-directed phase velocity and energy flux {dcu/dk > 0) and extraor¬ 
dinary, BW modes with contra-directed phase and group velocities {du/dk < 0). 
Some of them exhibit spatial dispersion which changes from positive to negative in 
different wavelength intervals whereas the others predominantly maintain negative 
dispersion. Frequencies and gaps between the modes can be tailored by changing 
shapes, sizes and spacing of the MM nano-building blocks so that phase velocities 
for the modes matching the photon energy conservation law would become equal 
as their group velocities remain contra-directed. Basically, different materials and 
geometries can be utilized to control spatial dispersion and losses in such MM with 
variable spatial dispersion. 


Second harmonic generation (SHC) is the basic fundamental coherent (i.e., 
phase-dependent) nonlinear optical propagation process which finds many prac¬ 
tically important application. Overall conversion efficiency is controlled by the 
intensity dependent local conversion rate and by the corresponding depletion of 
the fundamental beam. As shown in |TJ|2l|3l|ll|5] , such process exhibit unusual prop¬ 
erties in the BW settings which are significantly different from both BW optical 
parametrical amplification and SHC in ordinary materials. This paper is to ana¬ 
lyze such properties of SHC and fundamental differences between SHC in ordinary 
and backward-wave settings with focus on the pulse regime which is of particular 
interest. 
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2 Basic Equations 


As noted, the properties of SHG will be investigated for the cases where one of the 
coupled waves is ordinary and the other one is backward wave, in continuous-wave 
and pulse regimes, and compared with their ordinary counterparts. To achieve 
phase matching, phase velocities of the fundamental the SH waves must be co¬ 
directed. This means that in BW setting and pulse regime, the pulse of SH will 
propagate against the pulse of the fundamental radiation. Such extraordinary NLO 
coupling scheme dictates exotic properties of the frequency-doubling ultracompact 
NLO mirror under consideration, which are described below. One of our goals is to 
investigate the requirements for the magnitude of the metamaterial nonlinearity 
as well as for intensity of the control field(s) that would enable efficient functional 
properties of the indicated SHG mirror. Corresponding basic equations are as 
follows. 

Electric and magnetic components of the waves and corresponding nonlinear 
polarizations are defined as 


{£,71}j = Re{i5, H}j exp{i{kjZ — ojjt)}, 
{V, = Re {P, exp{i{kjZ — utjt)}, 

{P, ^ H}2{E, H}t, 

\P, Al f2 Xe,m^2\E, P fl, ^Xe,m,2 
Equations for amplitudes Ej can be written as 


dE2 1 dE2 . k20J2 , (2) ip2 r ■ A1 \ ^^2 Cl 

S 2 -Z -1-^ = -I - 2 E 1 exp (-lAkz) - -P-E 2 , 

oz V 2 at e 2 C^ ’ 2 

dEi 1 dEx -kiujio (2)* TP* TP f-Ai \ Cl 

si—-1- rr-=-’>■ - 2^'^Xei EiE2ex.p(iAkz) - —El. 

az vi at eic^ ’ 2 


( 1 ) 

( 2 ) 

(3) 

(4) 


(5) 

( 6 ) 


Here, Vi > 0 and 01^2 are group velocities and absorption indices at the correspond¬ 
ing frequencies, xi® = xi ^2 i® effective nonlinear susceptibility, Ak = k 2 — 2ki. 
Parameter Sj = 1 for ordinary, and Sj = —1 for backward wave. With ac¬ 
count for = n^{u)fc)'^, ni = 712 = S 2 y/e 2 ii 2 : we introduce amplitudes 

Cj = T^\ej\/kjEj, Oj = Cifeio, coupling parameters ae = •\/A:ifc 2 /|eie 2 |^ind 
g = seEio, loss and phase mismatch parameters ai ,2 = 01 , 2 ^ and Ak = Akl, slub 
thickness d = L/l, position ^ = zjl and time instant r = t/ At. It is assumed that 
EjQ = Ej{z = 0), I = viAr is the pump pulse length, Ar is duration of the input 
fundamental pulse. Quantities |aj|^ are proportional to the time dependent photon 
fluxes. Then Eqs. and are written as 


da2 vi da2 . 2 , ■ a 7 “2 

si^ + ^ = -i 2 c/*Zaia 2 exp - ^ai- (8) 

In the case of magnetic nonlinearity, xi^ = Xm^ 2 i equations for the magnetic 
components of the fields mj = Aj\pj\/kjHj, aj = mi/mio, take the form of (ff)) 
and ® where coupling parameters are = a/A:i fc 2 /|^i/r 2 1drrXg^ and g = eeHio. In 
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both cases, zq = g~^ is characteristic medium length required for signihcant NLO 
energy conversion. Only most favorable case of the exact phase matching Ak = 0 
and loss-free medium (ai = 02 = 0) will be considered below. 


3 Phase matched SHG of continuous wave in loss-free medium 

For the case of Z —>■ 00 , da^jdr —>• 0 and qi _2 = Ak = 0, Eqs. © and (|8]) reduce to 


sidai/dz = —iga\, s\da\ldz = —i2g*aia2- (9) 

Manley-Rove equation (photon conservation law) derived from Bps.© is 

S 2 d\a 2 f /dz -|- {\/2)s\d\a\\^/dz = 0. (10) 

3.1 SHG in ordinary NLO medium 

For the case of ordinary waves si = S 2 = 1 and 020 = 0 with account for aio = 1, 
one finds from Eq. m that 

2|a2|^ + |ai|^ = 1. (11) 

Solution to Eqs. m is found as [2113] 

2|u2|^ = tanh^ {V2gz), |ai|^ = sech^ {V2gz). (12) 


Eigures[T] (a) and (b) schematically show coupling geometry and energy fluxes 



Si 



§2 

^2 




Fig. 1 Coupling geometry and energy fluxes in ordinary, positive index, materials. 

S' 2,1 oc |a 2 ,i|^ across the ordinary, positive index material slab. According to 
Eqs. (HID and (11211 . sum of fluxes of photons huj 2 and of the pairs of photons 
huji is conserved across the medium. 
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Fig. 2 Coupling geometry and energy fluxes in negative-index metamaterials. 


3.2 SHG in the nanowaveguide with mixed, negative at ui and positive at 
2u) spatial dispersion 

Figures [2] (a) and (b) display coupling geometry and energy fluxes in the NLO 
slab in this case. In the case of backward fundamental and ordinary SH waves, 
Si = — 1 ,S 2 = 1, equations © and m dictates fundamentally different behavior. 
Equation dini) predicts 

\aif -2\a2f = \Bf, (13) 

where \B\'^ is a constant which, however, depends on the slab thickness and on 
strength of the input fundamental field. Equations © reduce to 

da^ldz = —igai, dai/dz = i2g*a\a2- (14) 

Besides the fact that the in this case the equations have different signs in the right 
sides, the boundary conditions for fundamental and SH waves must be applied to 
opposite edges of the slabs of thickness L: aio = 1, a 2 L = 0- Indicated differences 
give rise to fundamental changes in the solution to the equations for the amplitudes 
and, consequently, to the properties of the SHG as a whole. Solution to Eqs. m 
can be presented as P[3] 

V 2 a 2 = Bt&n[BV2g{L-z)], (15) 

ai = Bsec[BV2g{L-z)]. (16) 

As follows from (HU), quantity B presents transmitted fundamental wave, = B. 
It is found from Eq. (USD at 2 = 0 as: 

B = cos{BV2gL). (17) 

According to Eqs. m and (Hsi) - (dsi), difference between the numbers of pairs 
of photons hui\ and the number of photons hLL >2 is conserved so that the number of 
photons 1ul! 2 aX z = L is always equal to zero. The later requirement controls whole 
process of BWSHG. Such a dependence is in strict contrast with SHG in ordinary 
NLO materials, where the SH flux may exceed the fundamental one in some area 
inside the slab. In the case of BWSHG, its energy flux is directed against the 
fundamental one towards the areas with higher intensity of the fundamental wave 
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whereas remains lower everywhere inside the slab. As seen from Eq. m, the gap 
between the two curves, decreases and, hence, the conversion efficiency grows 
with gL = L/zo —> oo. Like in the preceding case, zq represent a characteristic 
medium length required for significant energy conversion to the SH. It shortens 
with growth of strength of the pump field. 



Fig. 3 Differences between SHG in ordinary and BW settings, (a). Energy ffuxes across the 
slab for fundamental (the descending dash-dotted blue line) and SH (the ascending dotted 
blue line), both are ordinary waves in a positive-index (PIM) material; and for backward-wave 
fundamental (the descending solid red line) and ordinary-wave SH (the descending dashed 
red line) in a negative-index (NIM) slab. gL = 1. (b).Output transmitted fundamental (the 
descending blue dash-dotted line) and SH (the ascending blue dotted line) sA, z = L, both are 
ordinary waves in a PIM; and transmitted backward-wave fundamental (the descending solid 
red line) flux at z = L and ordinary SH flux a,t z = 0 (the ascending dashed red line) in a 
NIM. 


3.3 Comparison of SHG for ordinary and backward fundamental waves 

Figure a) illustrate unparalleled properties of SH generation with backward 
waves as compared with its ordinary counterpart at similar other parameters 
for the particular example of pL = 1. As noted, the remarkable property is the 
fact that SH propagates against the fundamental wave and, therefore, metaslabs 
operates as a frequency doubling metamirror with reflectivity controlled by the 
fundamental wave. Further materials which support only ordinary waves and 
MM that support both ordinary and backward waves will be referred to as PIM 
and NIM. Figure EKb) compares output intensity of the SH and transmitted 
fundamental wave. Major important conclusion are as follows. It appears that the 
efficiency of SHG in ordinary settings exceeds that in the NIM metaslab at equal 
other parameters. At that, propagation properties of SH appear fundamentally 
different which holds promise of extraordinary applications. Intensity of SH is less 
than that of the fundamental wave across the NIM slab, whereas it can exceed 
that in the vicinity of exit from PIM slab [Fig. [Sja)] Quantum conversion grows 
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sharper with increase of intensity of fundamental beam and higher conversion at 
lower intensities occurs in PIM as compared with NIM. 

4 Comparison of ordinary and BW SHG in short-pulse regime 

In this section, comparative analysis of SHG in the ordinary and BW settings will 
be done for the case where intensity of input fundamental wave varies in time and, 
therefore, SHG is described by the system of partial differential equations © and 
©• The input pulse shape is chosen close to a rectangular form 



(18) 


where Sr is the duration of the pulse front and tail, and tq is the shift of the 
front relative to t = 0. Parameters Sr = 0.01 and tq = 0.1 have been selected for 
numerical simulations. Absorption is neglected (ai =02 =0). Phase velocities are 
equal Ak = 0. Group velocities are equal (V 12 ) and contra-directed. 

Unusual properties of SHG in NIMs in the pulsed regime stem from the fact 
that it occurs only inside the traveling pulse of fundamental radiation. Generation 
begins on its leading edge, grows towards its trailing edge, and then exits the fun¬ 
damental pulse with no further changes. Since the fundamental pulse propagates 
across the slab, the duration of the SH pulse may occur significantly longer than 
that of the fundamental one. Depletion of the fundamental radiation along the 
pulse and the overall conversion efficiency depend not only on maximum intensity 
of the input pulse, on matching of phase and group velocities of the fundamental 
and second harmonic, but on the ratio of the fundamental pulse length and slab 
thickness. Such properties are in strict contrast with that of SHG in PIM as il¬ 
lustrated in Figs. |4] and [5l Shape of the input fundamental pulse is given by the 
function Ti = |ai(r, 2 = 0)p/|aio|^ when its leading front enters the medium. The 
results of numerical simulations for the output fundamental pulse, when its tail 
reaches the slab boundary, is given by Ti = |ai(r, 2 = L)|^/|aio|^. Shape of the 
output pulse of SH, when its tail passes the slab’s edge at 2 = 0, are given by 
the function 772 = \a‘ 2 {'r,z = 0)|^/|aiop. Pulse energy are represented by the time 
integrated pulse areas Sj which varies across the slab. As seen from Figs. |4] (a) 
and (b), saturation is homogeneous across the pulse in PIM and shape of both 
fundamental and SH output pulses remain rectangular. On the contrary, shapes 
of the output fundamental and SH pulses in PIM are different and change with 
change of intensity of input fundamental pulse. Moreover, it appears that shapes of 
the output pulses varies with change of the input pulse length at unchanged other 
parameters as seen from Figs. [5] (a). Basic properties of pulse energy conversion 
across the NIM and PIM slabs qualitatively resemble those in continuous-wave 
regime [Figs. H] (c) and (d)]. Unparallel property in NIM is growth of pulse en¬ 
ergy conversion with shortening of pulse length at it constant instant intensity [cf. 
Figs, m (c) and [ 5 ] (b)]. Figures [5l[ a) and (b) correspond to the fundamental pulse 
four time shorter than the slab thickness. They show an increase of the conver¬ 
sion efficiency with increase of intensity of the input pulse. It is followed by the 
shortening of the SH pulse. 

Figures E] (c), (d) and [5] (b) satisfy to the conservation law in a loss-free 
metaslab: the number of annihilated pair of photons of fundamental radiation 
(Sio —S'ii)/2 is equal to the number of output SH photons 820 - Figures E] (a), (b) 
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(a) 


(b) 



Fig. 4 Comparison of pulse shapes and energy conversion at SHG in ordinary (PIM) and 
backward-wave (NIM) settings in a loss-free MM. The length of input pulse at the fundamental 
frequency is equal to the metaslab thickness. (a,b): Input rectangular Ti pulse shapes for the 
fundamental radiation; 772 - for SH. (c,d): Change of the pulse energy at the corresponding 
frequencies across the slab. Here, d = L/l, Si{z) and Sio = Si{z = 0) are fundamental pulse 
energy, 2 S 2 /S 10 is energy (photon) conversion efficiency per pulse. 


and [5Ka) prove that shapes and width of fundamental and generated SH pulses as 
well as the energy conversion efficiency to the reflected pulses at doubled frequency 
can be controlled by changing intensity and ratio of the the input pulse length to 
the metamaterial thickness (parameter d). 


5 Conclusion 

Comparative analysis of second harmonic generation in ordinary and backward- 
wave settings in continuous and pulse regimes is presented. Backward-wave regime 
is attributed to the nanomaterials which possess spatial dispersion of opposite signs 
at fundamental at second harmonic frequencies. This makes possible equal phase 
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Fig. 5 Effect of pulse width. Here input pulse duration is decreased in four times as compared 
with Fig. [4| (a) at the same peak intensity. Ti is pulse shape for the fundamental radiation, 
772 - for SH, Si{z) and ^lo = S\{z = 0) are fundamental pulse energy, 252/510 is energy 
(photon) conversion efficiency per pulse, d = L/l, 


velocities of the coupled waves whereas their energy fluxes appear contra-directed. 
Such deliberately engineered nanolayers can be viewed as nanowaveguides. It is 
shown that properties of second harmonic generation in ordinary and backward- 
wave settings are fundamentally different. In the latter case, metaslab serves as 
microscopic frequency doubling nonlinear-optical mirror which properties can be 
all-optically controlled. In pulse regime, energy conversion property are qualita¬ 
tively similar to those in continuous wave regime, whereas depend on the ratio 
of pulse length to the metaslab thickness. It is shown that width and shape of 
transmitted fundamental and generated second harmonic pulses can be controlled 
by changing intensity and width of input pulses. 
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